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INTRODUCTION 


1  . 

It  has  been  proposed  that  NO2  chemiluminescence  from  the  NO  +  0  and  NO 
*  O3  reactions  would  be  a  significant  source  of  atmospheric  background  IR 
radiance  in  a  disturbed  atmosphere . ^ ■ 2 ^  Predicting  this  radiance  requires 
a  reliable  model  of  N02  IR  emission  in  addition  to  knowledge  of  the  species 
concentrations.  Several  radiance  models  have  been  proposed,  each  of  which 
uses  different  assumptions  and  gives  significantly  different  results. 
Archer's  model^)  assumes  a  spectral  intensity  used  in  the  WOE  code,^3) 
based  on  "theoretical  considerations  of  energy  conservation."  This  spectrum 
gives  all  3  fundamental  bands  a  comparable  intensity,  and  also  includes 
numerous  overtone  and  combination  bands,  resulting  in  substantial  radiance 
throughout  the  LWIR.  A  second  model,  due  to  Kofsky  et  al.,^2,4^  is  based 
on  experimental  data  which  observes  roughly  a  10:1  ratio  of  v3  to  Vj  +  Vg 
intensity,  with  no  other  bands  included.  Finally,  a  more  sophisticated, 
room  temperature  model  for  v3  and  Vj+v3  chemiluminescence  has  recently  been 
developed  by  the  author. which  gives  excellent  agreement  with 
laboratory  data  at  pressures  in  the  0.02  -  2  torr  range. 

In  this  Report  a  simple  N02  IR  chemiluminescence  model  is  developed 
for  use  at  the  pressures  and  temperatures  found  in  the  upper  atmosphere. 
The  model  is  used  to  calculate  radiance  in  the  natural  atmosphere  in  the  80 
-  135  km  range,  where  comparison  can  be  made  with  the  predictions  of  Kofsky 
et  al.^2)  Since  the  radiance  is  proportional  to  the  NO  density,  the 
calculation  may  be  extrapolated  to  the  conditions  of  elevated  NO  found  in  a 
disturbed  atmosphere. 

The  current  model  is  based  on  the  author's  previous  work,^’®)  and 
gives  similar  results  at  laboratory  pressures.  At  very  low  pressures  there 
are  factor-of-3  differences  for  the  v3  and  Vj+v3  bands,  however,  and 
furthermore  the  current  model  includes  additional  bands  (v3  and  v2)  not 
included  previously. 
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2.  CURRENT  KNOWLEDGE  OF  THE  N02  CHEMILUMINESCENCE  MECHANISM 

As  discussed  elsewhere, <5*6)  available  data  on  N02  chemiluminescence 
is  consistent  with  the  following  mechanism: 

1.  The  NO  +  0  reaction  initially  yields  an  activated  complex  that 
appears  to  possess  energy  equipartition  between  its  electronic 
and  vibrational  modes.  This  complex  relaxes  either  radiatively 
(via  continuum  emission)  or  collisionally  to  the  ground 
electronic  state,  with  vibrational  equipartition  maintained  in 
either  case  as  a  consequence  of  chaotic  dynamics  in  the  excited 
electronic  (2B2)  state.  In  the  low  pressure  limit  the 
vibrational  distribution  following  continuum  emission  is 
characterized  by  an  average  energy  of  about  5000  cm-1.  These 
conclusions  are  derived  from  the  shapes  and  absolute 
intensities  of  the  continuum  and  IR  band  spectra  in  the  NO  +  0 
reaction. 

2.  The  NO  +  03  reaction  directly  yields  N02  in  vibrational 

equipartition.  The  nascent  state  distribution  is  close  to  a 
statistical  RRHO  distribution,  with  an  average  energy  of  about 
6000  cm-1.  These  conclusions  are  derived  from  trajectory 
calculations^ 7)  and  are  fully  consistent  with  the  shapes  and 
absolute  intensities  of  the  continuum  and  IR  band  spectra  in 
the  NO  +  03  reaction. 

3.  During  collisional  relaxation  of  vlbrationally  excited  N02, 
equipartition  between  the  Vj  and  v3  modes  is  approximately 
maintained,  resulting  in  a  -10:1  ratio  of  v3  to  Vj+v3  emission. 

This  indicates  that  Vj  and  v3  quanta  are  collisionally  removed 
with  comparable  rates.  This  behavior  is  consistent  with  either 
rapid  intermode  V-V  coupling,  or  similar  V-R.T  rates  for  Vj  and 
v3,  or,  as  suggested  by  studies  of  vibrational  deactivation  in 
other  triatomics , a  compromise  between  the  two 
mechanisms.  No  information  is  available  on  the  relaxation  of 
v2.  In  triatomics  such  as  03.^®^  C02,^9^  and  0CS,(]°)  the  rate 
of  intermode  coupling  is  roughly  comparable  to  the  rate  of 
relaxation  of  the  v2  mode,  leading  to  similar  effective 
quenching  rates  for  all  3  vibrational  modes.  Similar  behavior 
would  be  expected  for  N09 


4.  N02  reacts  with  atomic  oxygen  at  a  rate  of  0.9-1 .5xl0'ai  cm'3 
molec-1  sec-1,  depending  on  vibrational  excitation.  Using 
typical  N2,  02  and  0  atom  densities  in  the  100-135  km  altitude 
range,  where  N02  vibrational  emission  should  peak,  and  assuming 
N02  vibrational  quenching  rates  measured  in  the  laboratory  for 
N2  and  02,<11>12)  it  is  seen  that  most  of  the  collislonal 
removal  of  N02(v)  proceeds  via  reaction  with  0  rather  than 
physical  quenching.  This  means  that  even  if  the  assumption  o.,' 
equal  quenching  rates  for  the  3  vibrational  modes  is 
inaccurate,  the  Impact  on  the  column-integrated  atmospheric 
radiance  is  small. 

5.  The  Einstein  A  coefficients  for  emission  are  consistent  with 
harmonic  oscillator  scaling  rules;  i.e.,  A  -  v  for  fundamental 
bands.  In  this  approximation,  and  neglecting  overtone  and 
combination  bands,  the  rate  constant  for  radiative  removal  of 
energy  from  a  given  mode  is  simply  the  Einstein  coefficient  of 
the  v=l  fundamental. 


3.  MODEL  DESCRIPTION 


The  model  employs  equations  similar  to  Archer's  Eq.  (B-4),^1) 

I(j)  *  [NOHOHk2+2k3ifMi])*n)Qlj)  (1) 

1 

for  the  NO  +  0  reaction,  and 

I(j)  =  [N0][03]k14(j)Q(j)  (2) 

for  the  NO  ♦  0  reaction,  where  I(j)  is  the  radiance  of  band  j  in 
photons/cm3.  kj  is  the  NO  ♦  03  rate  constant,  k2  and  k3  are  the  2-body  and 
3-body  NO  +  0  rate  constants,  ♦  (j)  is  the  quantua  yield,  and  Q(j)  is  the 
quenching  factor,  given  by 

Q(J)  *  A  j  /  qd  [Mj  ]  +k0[0]  +A  j  )  (3) 

i 

where  Aj  is  the  Einstein  coefficient  of  band  j,  qj  is  the  vibrational 
quenching  rate  constant  for  species  Mj ,  and  kQ  is  the  rate  of  reaction  with 
0  atoas,  taken  as  lxlO-11  ca3  molec-1  sec-1. 

The  rate  constants  kj ,  k2  and  k3  have  roon  teaperature  values  given 
previously . The  temperature  dependences  are  taken  froa  Sharp^13)  for  the 
NO  +  0  reaction,  and  Borders  and  Birks^14^  for  the  NO  +  03  reaction. 

Assuming  that  5000  cm-1  of  vibrational  energy  is  deposited  in  N02  in 
both  reactions,  and  assuming  equiparti tion  among  the  3  vibrational  modes, 
the  quantum  yields  are  1.1,  2.2,  and  1.0  for  Vj,  v2,  and  v3,  respectively, 
qj  is  taken  as  5xl0-13  cm3  molec-1  sec-1  for  M  =  02  and  2xl0-13  for  M  = 
Ng .(11,12)  The  Einstein  coefficients  are  taken  as  130  sec-1  for  v3  (Ref. 
5)  and  0  02  sec-1  for  v2,  based  on  the  measured  band  strength .  ^ 13 ^  The 
Einstein  coefficient  for  the  very  weak  Vj  band  is  unknown,  since  a  reliable 


band  strength  is  not  available.  For  the  sake  of  argument  we  take  the  value 
0.01  sec'*,  which  probably  is  an  upper  limit. 

Unfortunately,  we  do  not  have  a  simple  way  of  estimating  effective 
quantum  yields  for  overtone  or  combination  bands.  For  the  present  we  shall 
take  the  Vj+v3  band  intensity  as  0.1  of  the  v3  intensity,  consistent  with 

-4  ./ 

both  experiments  and  modeling  at  laboratory  pressures.  This  ratio  might 
change  slightly  under  upper  atmospheric  conditions,  however.  We  assume 
that  emission  from  other  2-quantum  bands  is  much  smaller  than  v^Vg 
emission,  and  can  be  neglected. 

The  current  formulation  is  fully  consistent  with  the  previous  model 
calculations^5^  at  laboratory  pressures.  However,  it  yields  factor-of-3 
smaller  v3  and  Vj  +  v3  radiances  in  the  low  pressure  limit.  This  is  because 
the  previous  model  assumed  essentially  Instantaneous  Intermode  vibrational 
mixing,  causing  virtually  all  of  the  emitted  vibrational  energy  to  appear 
in  v3.  The  present  model  represents  the  opposite  point  of  view,  and 
assumes  that  the  energy  initially  apportioned  to  each  mode  (l.e.,  1/3  of 
the  available  vibrational  energy)  cannot  be  radiated  by  another  mode.  This 
also  appears  to  be  the  point  of  view  of  the  WOE  spectrum  used  by  Archer. 
Put  another  way.  the  previous  model  assumed  an  effective  v3  quantum  yield 
3x  as  large,  and  an  effective  Einstein  coefficient  1/3  as  large.  The 
neglect  of  lntermode  mixing  in  the  present  model  means  that  the  v3  radiance 
may  be  underestimated  and  that  the  predicted  Vj  and  v2  radiances  may  be 
overestimated  in  the  low  pressure,  low  0  atom  limit.  However,  as  explained 
in  paragraph  4.  of  Section  2,  the  column-integrated  N02  radiance  should  be 
insensitive  to  details  of  vibrational  relaxation. 


4 .  RESULTS 


Vj,  v2*  v3  an£*  Vj*v3  band  radiances  were  calculated  using  the 
altitudes  and  nighttime  0.  NO  and  03  densities  shown  in  Table  1 1  -  3  of 
Kofsky  et  al.^2)  The  N2  and  02  densities  and  temperature  were  taken  fro* 
the  U.  S.  Standard  Atmosphere . (*6 )  The  input  data  are  given  in  Table  1,  and 
the  computed  radiances  are  given  in  Table  2.  Integrated  zenith  and  95  km 
limb  radiances  were  computed  for  the  v3.  Vj+v3,  and  v2  bands,  and  are  given 
in  Table  3  in  Kofsky's  units  of  w/km2-sr.  For  the  v3  and  Vj+v3  bands  these 
values  turn  out  to  be  a  factor  of  3  to  5  lower  than  Kofsky's.  We  note  that 
Kofsky  assumed  slightly  different  v3  and  Vj+V3  quantum  yields,  and 
neglected  N02(v)  quenching  and  reactive  removal  as  well  as  temperature 
dependence  of  k2-  Both  studies  find  that  the  contribution  from  the  NO  +  03 
reaction  is  negligible  at  these  altitudes. 


TABLE  1.  ATMOSPHERIC  INPUT  DATA 


Altitude/km 

N?/cm"3 

02/cm-3 

0/cm  3 

0-,/cm"3 

NO/c*"3 

T/°K 

80 

3.0E14 

8.1E13 

3.E10 

2.E8 

8.E6 

199 

90 

5.5E13 

1  5E13 

2. Ell 

2.E8 

1  .E7 

187 

100 

9.2E12 

2.2E12 

5. Ell 

3.E7 

4  .  E7 

195 

110 

1  6E12 

2.6E11 

2. Ell 

2.E6 

5.E7 

240 

120 

3.7E11 

4 . 4E1 0 

1  .Ell 

2.E5 

3.E7 

360 

135 

8.7E10 

8.6E9 

4  .E10 

2.E4 

8.E6 

517 

to?. 


5.  DISCUSSION 


Although  the  present  Model  uses  equations  similar  to  Archer's,  and 
yields  a  similar  spectrum  In  the  collision-free  limit,  a  very  important 
difference  is  the  use  of  different  Einstein  coefficients  (hence,  different 
quenching  factors)  for  each  band.  This  results  in  spectra  which  are 
strongly  dominated  by  v3  and  Vj-v3,  as  predicted  by  Kofsky  et  al.  The 
WOE-type  spectrum,  which  shows  comparable  emission  in  all  bands,  is 
achieved  only  at  the  highest  altitudes,  where  the  intensity  is  weak  due  to 
the  falloff  of  NO  and  0  densities. 

A  qualitative  explanation  for  the  altitude  dependence  of  the 
vibrational  band  radiance  in  the  100-120  km  range  is  as  follows.  Because 
of  the  large  v3  Einstein  coefficient,  this  band  is  essentially  unquenched, 
hence  it  roughly  tracks  the  [0](N0]  product,  which  peaks  at  110  km.  On  the 
other  hand,  the  Vj  and  v2  bands  are  efficiently  removed  by  0  atoms,  hence 
these  bands  roughly  track  [NO],  which  is  nearly  constant  at  these 
altitudes.  At  lower  and  higher  altitudes  different  considerations  apply. 
At  135  km  all  the  bands  are  unquenched.  Below  100  km  all  the  bands  are 
quenched,  and  in  addition  the  3-body  contribution  to  the  NO  *  0  reaction  is 
important . 

Since  the  v3  and  Vj+v3  radiances  are  governed  by  parameters  which  are 
ultimately  tied  to  laboratory  data,  the  uncertainty  in  these  bands  is 
probably  comparable  to  the  factor  of  2  estimated  in  the  previous  model 
study  of  N02  vibrational  emission. ^  The  uncertainty  for  the  v2  band  is 
probably  larger,  since  this  band  has  not  been  observed  in 

chemiluminescence.  The  uncertainty  arises  mainly  in  our  estimate  of  the  v2 
quantum  yield,  which  may  lead  to  an  overall  uncertainty  of  perhaps  a  factor 
of  4.  the  lower  end  of  the  range  being  more  probable  than  the  higher  end 
The  Vj  band  radiance  is  even  more  uncertain  because  of  the  unknown  Einstein 
coefficient,  but  in  any  case  it  will  be  small  enough  to  be  of  negligible 
importance . 


The  v3  and  VJ  +  V3  band  radiances  nay  be  combined  with  experimental  or 
calculated  spectral  shapes  (e.g.,  Ref. 5)  to  yield  approximate  radiances  per 
■icron.  It  should  be  noted  that  these  spectra,  particularly  V2+V3,  are 
broadened  and  red  shifted  from  the  corresponding  cold  bands.  The  V£  band 
should  exhibit  similar  effects. 
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